INFINITE BAR-JOINT FRAMEWORKS, CRYSTALS 
AND OPERATOR THEORY 



J.C. OWEN AND S.C. POWER 

Abstract. A theory of flexibility and rigidity is developed for 
general infinite bar-joint frameworks (G,p). Determinations of 
nondcformability through vanishing flexibility are obtained as well 
as sufficient conditions for deformability. Forms of infinitesimal 
flexibility arc defined in terms of the operator theory of the as- 
sociated infinite rigidity matrix R(G,p). The matricial symbol 
function of an abstract crystal framework is introduced, being the 
matrix- valued function on the (i-torus representing R(G,p) as a 
Hilbcrt space operator. The symbol function is related to infini- 
tesimal flexibility, deformability and isostaticity. Various generic 
abstract crystal frameworks which are in Maxwellian equilibrium, 
such as certain 4-regular planar frameworks, are proven to be 
squarc-summably infinitesimally rigid as well as smoothly deformable 
in infinitely many ways. The symbol function of a three-dimensional 
crystal framework determines the infinitesimal wave flexes in mod- 
els for the low energy vibrational modes (RUMs) in material crys- 
tals. For crystal frameworks with inversion symmetry it is shown 
that the RUMS appear in surfaces, generalising a result of F. Weg- 
ner [32] for tetrahedral crystals. 
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1. Introduction 



Infinite bar-joint frameworks appear frequently as idealised models 
in the analysis of deformations and vibration modes of amorphous and 
crystalline materials. See [11], [18], [6], [H], [32] and [36] for example 
and the comments below. Despite these connections there has been no 
extended mathematical analysis of such models. Notions of rigidity, 
flexibility, deformability, constrainedness, independence and isostatic- 
ity, for example, are usually employed either in the sense of their usage 
for a finite approximating framework or in a manner drawn from ex- 
perience and empirical fact in the light of the application at hand. It 
seems that a deeper understanding of the models is of considerable in- 
terest in its own right and that a mathematical development may prove 
useful in certain applications. In what follows we shall provide formal 
definitions of the terms above in quite a wide variety of forms and we 
examine some of their inter-relationships and manifestations. 

Suppose that one starts with a flexible square bar-joint framework 
in two dimensions and that this is then extended periodically to create 
an infinite periodic bar-joint network. Is the resulting assemblage, with 
inextendible bars, continuously flexible in two dimensional space? A 
moment's reflection reveals a proliferation of flexibility, such as sheering 
motions with one half of the network fixed. 

However such movement is dramatically infinite and a natural sec- 
ond question is whether for periodic frameworks such as these there 
are flexes for which the total joint movement is finite. The less obvious 
negative answer in these cases offers some satisfaction in that it is con- 
sistent with so-called Maxwell counting in the case where the average 
number of degrees of freedom of the joints matches the average number 
of constraints per joint. (See Theorem 5.2 and Corollary 5.3.) 

On the other hand if, for the grid example, one rigidities alternate 
squares by adding diagonal bars, as shown, then the resulting structure 
of corner linked rigid squares remains properly flexible, although now 
uniquely so. In fact the unique flex has an affinely contracting character 
(See Definition 4.3) with alternating rotation of the squares. 
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More generally the flexibility of polytope networks in two and three 
dimensions continues to be of interest in the modeling of crystals and 
amorphous materials, especially with regard to their low frequency vi- 
brational modes. Such modes appear, for example, in higher order 
symmetry phases of tetrahedral crystals and are referred to as rigid 
unit modes (RUMs). Indeed in the paper of Giddy et al [H] the alter- 
nating flex of the squares framework above has been associated with 
vibrational modes in perovskite. See also Hammond et al [18], Wegner 
[32], as well as Goodwin et al [H] for a useful overview. At the same 
time, in the modeling of amorphous materials, such as glasses, there is 
interest in understanding the critical probabilities that guarantee flexi- 
bility and rigidity for classes of randomly constructed frameworks. See, 
for example, Chubynsky and Thorpe [B] for the recent determination 
of such probabilities in simulation experiments. 




Figure 1. The grid framework in the plane, Q%2. 




Figure 2. The corner-joined squares framework, Q, 



Formally, a framework in R d (or bar-joint framework, or distance- 
constraint framework) is a pair Q = {G,p) where G = (V, E) is a simple 
connected graph and p = (pi,P2, ■ ■ •) is a framework vector made up of 
framework points pi in M. d associated with the vertices v±, V2, ■ ■ ■ of V. 
The framework edges are the (closed) line segments [pi,Pj] associated 
with the edges E of the graph G = (V, E). As the ellipsis suggest, we 
allow G to be a countable graph. We shall also define a crystal frame- 
work C as a framework with translational symmetry which is generated 
by a connected finite motif of edges and vertices. (See Definition 4.2.) 
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When G is finite and the framework points are generically located 
in M 2 then a celebrated theorem of Laman [2TJ, well-known in struc- 
tural engineering and in the discrete mathematics of rigidity matroids 
[To] , gives a simple combinatorial criterion for the minimal infinitesimal 
rigidity of the framework; the graph itself satisfies Maxwell's counting 
rule 2\V\ — \E\ = 3, and subgraphs G' = (V',E') must comply with 
2\V'\ — \E'\ > 3. This is a beautiful result since the rigidity here is the 
noncombinatorial requirement that the kernel of an associated rigid- 
ity matrix R(G,p) has the smallest dimension (namely three) for some 
(and hence all) generic framework. On the other hand frameworks 
with global symmetries, or even with "symmetric elements" (such as 
parallel edges) are not generic, that is, algebraic dependencies do exist 
between the framework point coordinates. Such frameworks arise in 
classical crystallography on the one hand and in mathematical models 
in structural engineering and in materials science on the other. See, for 
example, Donev and Torquato [7], Hutchinson and Fleck [19], Guest 
and Hutchinson [16] and various papers in the conference proceedings 

The present paper develops two themes. The first concerns a math- 
ematical theory of deformability and rigidity for general infinite frame- 
works, with frequent attention to the case of periodic frameworks. 
There is, unsurprisingly, a great diversity of infinite framework flex- 
ing phenomena and we introduce strict terminology and some methods 
from functional analysis to capture some of this. In the second theme 
we propose an operator theory perspective for the infinitesimal (first 
order) flexibility of infinite frameworks. 

Particularly interesting classes of infinite frameworks, from the point 
of view of flexibility, are those in the plane whose graphs are 4-regular 
and those in three dimensions whose graphs are 6-regular. In this 
case the graphs are in Maxwell equilibrium, so to speak, and so in a 
generic framework realisation any flex must activate countably many 
vertices. This is also the case for various periodic realisations such 
as the kagome framework, Qkag, formed by corner-joined triangles in 
regular hexagonal arrangement, and frameworks in three dimensions 
formed by pairwise corner-joined tetrahedra. Despite being internally 
rigid in this way (Definition 2.18 (vi)) these frameworks admit diverse 
deformations. For example we note that the kagome framework admits 
uncountably many distinct deformations and in Theorem 4.4 we note 
that Z d -periodic cell-generic grid frameworks in M. d admit deformations 
associated with affine transformations. 

A significant phenomenon in the infinite setting is the appearance of 
vanishing flexibility. This means, roughly speaking, that the framework 
is a union of finite flexible subframeworks but the extent of flexibility 
diminishes to zero as the size of these subframeworks increases, so that 
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the infinite assemblage is inflexible. Elementary examples were indi- 
cated in [25] but we give more subtle examples here which are due to 
flex amplification at second order distances through concatenation ef- 
fects. In particular there are bounded infinitesimal flexes in periodic 
frameworks that admit no continuous extensions and which do not arise 
as the derivative of a smooth deformation. We also note that there are 
Z 2 -periodic crystal frameworks which are somewhat paradoxical, being 
indeformable despite the flexibility of all supercell subframe works. On 
the other hand, in the positive direction, in Theorem 12.201 we give a 
general result which identifies a uniform principle for the existence of a 
deformation. The proof uses the Ascoli-Arzela theorem on the precom- 
pactness of equicontinuous families of local flexes. It remains an inter- 
esting open problem to determine necessary and sufficient conditions 
for the rigidity and bounded rigidity of periodic planar frameworks. 

The operator theory perspective for frameworks was suggested in 
[25] as an approach to a wider understanding of infinitesimal flexibility 
and rigidity. In this consideration the rigidity matrix is infinite and de- 
termines operators between various normed sequence spaces associated 
with nodes and with edges. Also, in [27] we have given a direct proof 
of the Fowler-Guest formula [2j for symmetric finite frameworks which 
is based on the commutation properties of the rigidity matrix as a lin- 
ear transformation and this adapts readily to the infinite case and the 
rigidity operators of crystal frameworks. Indeed, translational symme- 
try ensures that the rigidity matrix R(G,p) intertwines the coordinate 
shift operations. We consider square summable flexes and stresses and 
for distance regular locally finite frameworks R(G,p) is interpreted as a 
bounded operator between Hilbert spaces. Also, enlarging to complex 
Hilbert spaces the Fourier transform J r R(G,p)J r_1 is identified as a 
multiplication operator 

M$ : L 2 (T d ) <g> C n -> L 2 (Y d ) <g> C m 

given by an m x n matrix-valued function <&(z) on the rf-torus. The 
function $ for C is referred to as the matricial symbol function as- 
sociated with the particular generating motif. The terminology and 
notation is borrowed from standard usage for Toeplitz operators and 
multiplication operators (see [3] for example). Many aspects of in- 
finitesimal flexibility and isostaticity are expressible and analysable 
in terms of the matricial symbol function and its associated opera- 
tor theory. For example a straightforward consequence of the operator 
theoretic approach is the square summable isostaticity of various non- 
degenerate regular frameworks that satisfy Maxwell counting, such as 
grid frameworks and the kagome framework. 

An explicit motif-to-matrix function algorithm is given for the pro- 
gression 

C = (G,p) ->■ R(G,p) -> $(z). 
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Furthermore the identification of infinitesimal periodic-modulo-phase 
flexes and their multiplicities is determined by the degeneracies of 3>(z) 
as z ranges on the <i-torus. In particular, the function 

H(z) : = dimker$(» : T d -> Z. 

gives a determination of the mode multiplicity of periodic-modulo- 
phase infinitesimal flexes. 

In the discussions below we are mainly concerned with properties 
of mathematical bar-joint frameworks. (The framework bars are in- 
destructibly inextensible, the joints are located deterministically, they 
maintain perfect, frictionless fit and may even coincide.) Nevertheless, 
analysis of matricial symbol functions and their degeneracies are par- 
ticularly relevant to the description and analysis of Rigid Unit Modes 
in material crystals. We show that for crystal frameworks with inver- 
sion symmetry the set of RUMS is a union of surfaces. This generalises 
and provides an alternative perspective for a recent result of Wegner 
[32] for tetrahedral crystals. 

Operator theory methods have proven beneficial in many areas of 
mathematics and applications. This is evidently the case for multi- 
linear systems theory and in control theory for example. Infinite rigid- 
ity matrix analysis seems to possess some similitudes with these areas 
and it seems to us that here too the operator turn will be a useful one. 

The development is as follows. Section 2 gives a self-contained ac- 
count of continuous flexibility and rigidity, and vanishing flexibility and 
one-sided flexibility is proven for various periodic infinite frameworks. 
Forms of flexibility, such as bounded flexes, square summable flexes, 
summable flexes and vanishing flexes are defined and determined for 
some specific examples. Sufficient conditions are obtained for the ex- 
istence of a smooth flex and a flex extension problem for generic finite 
frameworks is posed. A positive resolution of this problem would pro- 
vide a natural extension of Laman's theorem to infinite frameworks. 
In Section 3 we consider infinitesimal theory for general infinite frame- 
works and determine a number of rigidity operators and their flex and 
stress spaces. (The topic is taken up in more detail for crystal frame- 
works in Section 5.) In Section 4 we consider (abstract) crystal frame- 
works in two or three dimensions. These are generated by a motif and 
a discrete translation group. Various forms of deformations are con- 
sidered, such as strict periodic flexibility, flow-periodic flexibility and 
flexes with reduced periodicity and symmetry. Also we indicate the 
flat torus model for crystal frameworks and recent results of Ross and 
Whiteley in this direction, including a periodic analogue of Laman's 
theorem. In the final section we consider the matricial symbol function 
approach and various examples. In particular we determine the (unit 
cell) infinitesimal wave flex multiplicities for the kagome net framework 
by factoring the determinant of the matricial symbol function. 
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2. Infinite Bar-joint Frameworks 

In this section we give a self-contained rigourous development of in- 
finite frameworks and examine the nonlinear aspects of their flexibility 
by continuous deformations and their associated rigidity. In the next 
section we consider infinitesimal flexibility and rigidity in a variety of 
forms. 

2.1. Continuous flexibility and rigidity. We first define continuous 
flexes and continuous rigidity. The latter means, roughly speaking, that 
the framework admits no proper deformations that preserve the edge 
lengths. The definition below gives straightforward generalisations of 
terms used for finite frameworks. In that case we note that a continuous 
flex is often referred to as a finite flex while in engineering models it is 
referred to finite mechanism. 

Unless we specify otherwise we shall assume that the frameworks 
under consideration are proper in that the framework points do not lie 
on a hyperplane in the ambient space M. d and that the framework edges 
[Pi,Pj] have nonzero lengths |pj — Pj\. 

Definition 2.1. Let (G,p) be an infinite framework in IR 2 , with con- 
nected abstract graph G = (V, E), V = {v\, v 2 , . . . } andp = (pi,£>2> • • • )• 

(a) A base- fixed continuous flex, or, simply, a flex of (G,p), is a func- 
tion 

p{t) = (pi(t),p 2 (*), • • • ) from [0, 1] to Ylv ^ with the following prop- 
erties; 

(i) p(0) = p, 

(ii) each coordinate function Pi : [0, 1] — > M. 2 is continuous, 

(iii) for some base edge (v a ,v b ) with \p a — p b \ ^ 0, p a (t) = p a (0) and 
Pb(t) = Pft(0) for all t, 

(iv) each edge distance is conserved; \pi(t) — Pj(t)\ = \pi(0) — Pj(0)\ 
for all edges (vi,Vj), and all t, and 

(v) p{t) 7^ p for some t G (0, 1]. 

(b) The framework (G,p) is flexible, or more precisely, continuously 
flexible, if it possesses a base-fixed continuous flex. 

(c) The framework (G,p) is rigid, or continuously rigid, if it is not 
flexible. 
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Similarly one defines base-fixed continuous flexes and continuous 
rigidity for proper frameworks in M. d by replacing a base edge by an 
appropriate set of framework points with maximal affine span. 

The simplest kind of continuously rigid framework in the plane is 
one which is a union of continuously rigid finite frameworks. In partic- 
ular the following theorem follows simply from the theorem of Laman 
indicated in the introduction. 

Theorem 2.2. Let G be a connected graph which is the union of a 
sequence of finite Laman graphs. Then every generic realisation (G,p) 
in the plane is continuously rigid. 

By generic, or, more precisely, algebraically generic, we mean, as 
is usual, that the coordinates of any finite set of framework points is 
algebraically independent over the rational numbers. Unlike the case 
of finite frameworks it is possible to construct two generic frameworks 
with the same abstract graph one of which is flexible and one of which 
is rigid. Accordingly it seems appropriate to formulate the following 
definition to extend the usual usage. 

Definition 2.3. An infinite simple connected graph G is said to be 
rigid, or generically rigid, for two dimensions, if every generic frame- 
work (G,p) in the plane is rigid. 

Note that if G is rigid and if if is a containing graph for G with 
the same vertex set then every generic framework (H,p) in the plane 
is rigid. 

It seems likely that the converse to the theorem above holds. That is, 
if H does not contain a sequentially Laman graph (in the sense below) 
with the same vertex set, then H is not generically rigid. We comment 
more on this later in Section 2.6. 

Rigidity and flexibility are properties of the entire framework and it 
is such entire features and their inter-relationships that are of primary 
interest in what follows. One would like to understand the relationship 
with small scale or local structure, such as local counting conditions 
and local connectivity. Additionally, as above, one would like to relate 
entire properties to sequential features that pertain to an exhausting 
chain of finite subframeworks and for this the following definition is 
helpful. 

Definition 2.4. If P is a property for a class of finite, simple, connected 
graphs (resp. frameworks) then a graph G (resp. framework Q = 
(G,p)) is sequentially P or a-P if G is the union of graphs in some 
increasing sequence of vertex induced finite subgraphs G\ C G% C . . . , 
and each graph Gk (resp. framework (Gk,p)) has property P. 

For example, we may refer to an infinite graph as being a-Laman, or 
a-Laman— 1 and an infinite framework as being cr-rigid. To say that an 
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infinite framework is a-flexible, or sequentially flexible, is rather vac- 
uous since it usually prevails for trivial reasons. (One can construct 
countably infinite edge complete frameworks that may fail to be so but 
all our examples are sequentially flexible.) The more interesting prop- 
erty is the failure of sequentially rigidity as in the following definition. 

Definition 2.5. A framework Q is said to have vanishing flexibility if 
it is continuously rigid but not a-rigid. 

An important topic in finite rigidity frameworks is that of global 
rigidity, also called unique rigidity. This holds, in the two dimensional 
setting, if there is, up to congruency, only one embedding of the frame- 
work in R 2 . This implies, for example, that these frameworks admit 
no foldings, and indeed have unique diagrams up to rotations and re- 
flection. One can extend the term to infinite frameworks but we do 
not consider this issue here at all. One might be tempted to say that 
a rigid infinite framework, especially one with vanishing flexibility, is 
globally rigid, but we refrain from doing so because of conflict with this 
usage. 

The term "global" for global rigidity is natural since rigidity for 
finite frameworks is equivalent to the "local" property that "nearby 
equivalent frameworks are congruent". That is, if there exists e > 
such that if (G,p') is a finite framework with \pi — p'j\ < e, for all i, and 
if {G,p') is equivalent to (G,p), in the sense that corresponding edges 
have the same length, then (G,p') and (G,p) are congruent. See, for 
example, Gluck [12] and Asimow and Roth [TJ. 

An infinite simple graph G is locally finite if for every vertex v there 
are finitely many incident edges. Amongst such graphs are those for 
which there is an upper bound to the degree of the vertices, as in the 
case of the graphs of crystal frameworks. Within this class a graph G 
is said to be r -regular if every vertex has degree r. We remark that the 
theory of tilings provides a wealth of examples of planar frameworks 
which are 4-regular. 

Remark 2.6. In what follows we consider only locally finite frame- 
works. Without this assumption it is possible to construct quite wildly 
flexing planar linkages. In fact, given a continuous function / : [0, 1] — > 
IR 2 one can construct an infinite linkage, in the sense of the definition 
below, and a base-fixed flex p(t) with a motion p v (t) for a particular 
vertex v that is equal to /(£). This includes the possibility of space 
filling curves. This is a consequence of a continuous analogue of a well 
known theorem of Kempe which asserts that any finite algebraic curve 
in the plane can be simulated by a finite linkage. For more details see 
Owen and Power [26]. 
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2.2. Linkages. The removal of a framework edge from a rigid frame- 
work may result in flexibility which is, roughly speaking, of a one- 
dimensional nature. We reserve the term linkage for such a mathe- 
matical object, which we formally specify in the next definition. We 
remark that finite frameworks are also referred to as linkages, particu- 
larly when they are flexible, perhaps with several degrees of freedom, 
but this should not cause confusion. 

A two-sided continuous flex p{t) of (G, p) is defined as above but for 
the replacement of [0, 1] by [—1, 1]. The following formal definition of 
an infinite linkage reflects the fact that the initial motion of a base- fixed 
linkage is uniquely determined by the angle change at any flexible joint. 

Definition 2.7. A linkage in M? is a finite or infinite connected frame- 
work Q = (G,p) in M 2 for which there exists a continuous two-sided 
base- fixed flex p(t) with framework edges [pi,Pj], [Pj,Pk] such that the 
cosine angle function 

g(t) = (pi{t)- Pj (t),p k {t)- Pj (t)) 

is strictly increasing on [—1,1], and such that p(t) is the unique two- 
sided flex q(t) of Q with qi(t) = pi(t), for I = i, j, k. 

Many interesting finite linkages were considered in the nineteenth 
century in connection with mechanical linkages. See, for example, 
Kempe [20]. Note however that the definition is liberal in allowing 
coincident joints and self-intersecting flexes. Also the definition refers 
to local deformation behaviour and this does not rule out the possi- 
bility of bifurcations occurring in a parameter extension of the given 
flex. 

It is a simple matter to construct diverse infinite linkages by tower 
constructions or progressive assembly. (See, for example, the Cantor 
tree frameworks of [22].) However, some such constructions lead to 
frameworks with vanishing flexibility and so are not linkages in this 
case. An elementary illustration is given in Figure 4 wherein a two-way 
infinite rectangular strip linkage is augmented by adding flex-restricting 
cross braces in an alternating fashion. If the brace lengths tend to 
the diagonal length from above then the infinite framework is rigid. 
Evidently in this case the triangle inequality is playing a role in isolating 
one real solution to the solution set V(G,p) defined below. One can 
also construct examples in which this isolation is less evident, with all 
joint angles bounded away from zero and ir/2 for example. 

A more interesting and subtle form of vanishing flexibility is due to 
progressive flex amplification rather than local flex restrictions. Roughly 
speaking, if a small flex is initiated at a particular joint and the flex 
propagates in some amplifying manner, then the triangle inequality at 
some far remove may prohibit any further increase. If the framework 
is infinite then no local joint flex may be possible at all. The strip 
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framework of concatenated levers in Figure 5 gives an example where 
the amplification is evident, while the rigidity of the strip framework 
in Figure 6 and the trapezium strip in Figure 7 is less evident. The 
lever framework has a natural infintesimal flex, in the sense of Section 
3, which is unbounded. The corresponding flexes of Figures 6 and 7 
however, are bounded with amplification unfolding as a second order 
effect. This is proven in the next subsection. 




Figure 3. An infinite rectangle strip linkage. 




Figure 4. A restricted rectangle strip. 




Figure 5. Rigid but not cr-rigid. 




Figure 6. Rigid but not cr-rigid. 



It is a straightforward matter to incorporate the vanishing flexibility 
of the strips above as subframeworks of a Z 2 -periodic framework. This 
process is indicated in Figures 9, 10 and 11 below, where the infinite 
frameworks are determined as the periodic extensions of the given unit 
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Figure 7. Rigid but not a-rigid. 




Figure 8. A periodic half-strip which is only right-flexible. 

cell. Figure 9 shows a linkage formed as a "fence lattice" composed 
of infinite horizontal and vertical a-rigid bands. Figure 10 shows an 
analogue where the infinite bands have been replaced by rigid strip 
frameworks. Figure 11 is an elaboration of this in which cross braces 
have been introduced to remove the flexibility. (Only one edge is needed 
for this whereas the example given is periodic.) The additional degree 
2 vertex in the cell ensures that the framework is not cx-rigid, while the 
infinite bands remain vanishingly rigid. 

Let us note that for the framework in Figure 11, with its curious 
mixture of rigidity and flexibility, one can add any finite number of 
additional degree 2 vertices without changing the rigidity of the frame- 
work. In particular we have a construction that proves the following 
proposition. 




Figure 9. Unit cell for a "fence lattice" linkage. 



Proposition 2.8. Let c > 2. Then there is a 1? -periodic framework in 
M? which is rigid, which is not o-rigid and for which the average vertex 
degree is less than c. 
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Figure 10. Unit cell for a modified fence lattice linkage. 



Figure 11. Unit cell for a rigid periodic framework 
which is not a-rigid. 

One can readily extend this fanciful idea in various ways to obtain 
such structures in higher dimensions. For example, start with a one- 
dimensionally periodic a-rigid girder in 3D and augment it with trapez- 
ium "tents" of alternating height to creates vanishingly rigid girders. 
Also periodically interpolate any number of degree two vertices into the 
tent top edges without removing the vanishing flexibility. Join infin- 
itely many such component girders periodically at appropriate tent-top 
edges to create a fence framework and add linear jointed cross braces 
to create, finally, a 2D periodic grid which is continuously rigid in 3D, 
which is not er-rigid and which has average coordination number arbi- 
trarily close to two. 

2.3. Relative rigidity and the extension of flexes. For a finite or 
infinite framework Q = (G,p) in M. 2 define the function 

f G : n R2 ri M ' fM = (I® - ^i 2 )e= ( ^)- 

V E 

This is the usual edge function of the framework and depends only on 
the abstract graph G. 

Definition 2.9. The solution set, or configuration space, of a frame- 
work Q = (G,p), denoted V(G,p), is the set /g 1 (/g(p))- This is the set 
of all framework vectors q for G that satisfy the distance constraints 
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equations 

\li ~ Qj? = \Pi ~Pj\ 2 i for a U edges e = (v h Vj). 

In general the solution set of an infinite framework need not be a real 
algebraic variety even when it is "finitely parametrised". In less wild 
situations it can be useful to relate V(G,p) to the algebraic variety 
V (H , it h (p)) associated with a finite subgraph H of G, or with an 
elementary subgraph such as a tree, or even a set of vertices. 

Definition 2.10. An infinite bar-joint framework (G,p) in M. d is rigid 
over a subframework (H, tth{p)) if every continuous flex of (G,p) which 
is constant valued on (H,tth{p)) is constant. Similarly, if H is a sub- 
graph of a countable connected simple graph G then G is rigid over 
H, or generically rigid over H if, for every generic frameworks (G,p), 
every continuous flex of {G,p) which is constant- valued on (H,ith(p)) 
is constant. 

We may also form the following associated notions. 

Definition 2.11. An infinite framework (G,p) in M. d is finitely deter- 
mined if it is rigid over (H, p) for some finite subgraph H and is finitely 
flexible if it is flexible and finitely determined. 

Finite flexibility in the sense above is a strong property in which 
paths from p in the solution set V(G,p) are determined near p by 
the finite algebraic variety V(H,p). Note that the term "infinitely 
flexible" is not appropriate to describe a flexible framework which is 
not finitely flexible since it is possible to construct linkages, in our 
formal sense, which are not finitely flexible. This is the case for the 
periodic framework in Figure 12 which is a linkage because of partial 
vanishing flexibility. 



Figure 12. A linkage which is not finitely determined. 



An important class of frameworks which appear in mathematical 
models are those that are distance regular. 

Definition 2.12. A framework Q = (G,p) in M. d is distance-regular if 
there exist < m < M such that for all edges 

m < \Pi — Pj I < M. 
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For such a framework (G,p) it is natural to consider the nearby 
frameworks with the same graphs but with slightly perturbed frame- 
work points (and therefore edge lengths). If a property holds for all 
such frameworks, for some perurbation distance e then we call such a 
property a stable property for the the framework. 

Formally, an e-perturbation of a distance regular framework Q = 
(G,p) is a framework Q' = (G,p') for which \p% — p'i\ < e, for all cor- 
responding framework points. Recall that a finite framework in M. d is 
said to be e-rigid if it is congruent to every equivalent e-perturbation. 
Let us say that a general framework is perturbationally rigid if it is 
e-rigid for some e. It is a well-known fact that perturbational rigidity 
and rigidity are equivalent in the case of algebraically generic finite 
frameworks [T], [12]. However, it is straightforward to see that this 
equivalence thoroughly fails for general infinite frameworks (see [25]). 

Definition 2.13. Let Q be a distance-regular framework. Then Q is 
stably rigid (resp stably flexible) if it is rigid (resp. flexible) and for 
sufficiently small e > every e-perturbation of Q is rigid (resp. flexible). 

Likewise, if P is any particular property of a distance-regular infinite 
framework then we may say that Q is stably P if, for some e > 0, the 
property P holds for all e-perturbations. 

Proposition 2.14. The periodic trapezium strip frameworks, with al- 
ternating unequal heights a > b > o, are rigid. In particular the rec- 
tangle strip linkage (of Figure 7) is not stably flexible. 




Figure 13. 



Proof. Figure 13 shows the displacement of a double trapezium to the 
right. Let the three vertical bar lengths be a,b and a units with a > 
b > 0. The displaced position has angles A,B,C at the base line 
and angles D, E, F, G occurring relative to the trajectory tangents of 
the displaced vertices. For a subsequent incremental change 8 A, with 
resulting incremental changes 8B, 8C, 8D, 8E, 8F, 8G we can see from 
simple geometry that to first order 

a(6 A cos D) = b(8B cos E). 
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Suppose now that A is regarded as a specialisation of the input angle 
a with resulting output angle (3 = (3(a), so that at a = A we have 
(5(a) = B. Then 

dj3 , SB a cos D 

da sa^o 5A bcosE' 

Similarly, with angle C regarded as the output angle j(B) for the angle 
transmission function 7 = j(/3), we have 

g?7 b cos F 

d/3 a cos G 

and so 

d 1 dj. d(3 J cos -D cos F 

da a ~ A d[5 da a ~ A cos E cos G 
Note that since B > A we have also D > E, and since B > C we have 
F > G, from which it follows that both ratios above are less than one. 
Thus certainly < 7' (a) < 1 for < a < a% where a\ is the first 
positive angle for which j'(ai) = 0. 

It follows, from the mean value theorem, that the double trapezium 
angle transmission function is an increasing differentiable function with 

7(0) = 0, < 7(a) < a, for < a < an. 

It follows immediately that the right-semi-infinite trapezium strip is 
right flexible. 

We let A = 7(011), which we refer to as the locking angle. Note the 
second trapezium of the double admits no increase of this angle. In 
view of the above we have A < a\. 

Suppose now that Q is the two-way-infinite trapezium strip, with 
a ^ b. Let p(t) be a flex and suppose that for a fixed framework edge 
with length a the angle A is greater than zero for some time ti > 
and that t± is the first such time. Then certainly < A < A. Note 
that A_ n = 7~ n (A), n = 1,2,... are the angles of the edges of length a, 
counted off to the left. In view of the function dominance < 7(a) < a 
it follows that A_ n > a\ for some n, which is a contradiction. □ 

The argument above also shows that the semi-infinite trapezium strip 
framework of Figure 7 has a continuous flex but no two-sided flex. 

Remark 2.15. It seems to be of interest to analyse strip frameworks in 
further detail. For example, a trapezium strip framework is not stably 
rigid, despite the apparent "robustness" of the argument above. To see 
this use surgery in the following way. Remove one cross bar, then push 
the rightmost semifinite strip to the right, by an angle perturbation 
A = e > 0. Now insert a replacement bar of the required length. One 
can flex the resulting structure towards the left to restore the position 
of the right hand strip. Indeed, this is all the flexibility the framework 
has. The possible flex of an e perturbation, such as the one described, 
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seems to be of order e and so there does seem to be "approximate 
rigidity" . 

Remark 2.16. Consider a periodic trapezium grid framework Qtrap 
obtained by perturbing Q%2 by adding a fixed small positive value to 
the y coordinate of the framework points pij for the odd values of i 
and j. It can be shown that this framework is rigid over any linear 
subframework. This contrasts with the grid framework itself which is 
freely flexible over its x and y axes in the following specific local sense: 
all sufficiently small flexes of the subframework extends to a flex of the 
whole framework. On the other hand note that Theorem 4.4 shows 
that Qtrap is deformable. 

2.4. Forms of flexibility. It seems to be a fundamental and interest- 
ing issue to determine the ways in which infinite bar-joint frameworks 
are rigid or continuously flexible. In this section we give some further 
definitions, we give sufficient conditions for the existence of a proper flex 
and we contemplate a plausible infinite framework version of Laman's 
theorem. 

Flexes are often infinitely differentiable or smooth in the sense of 
the following formal definition. This is the case for example, for the 
"alternation" flexes of Q sq and Qk ag - 

Definition 2.17. A continuous base-fixed two-sided flex p(t) : t 6 
[—1,1] of a framework (G,p) in M. d is a smooth flex if each coordinate 
function Pi(t) is infinitely differentiable. 

The smoothness of a flex is a local requirement whereas the following 
terms impose various increasing forms of global constraint. In partic- 
ular rotational flexes of infinite frameworks with unbounded diameter 
are not bounded flexes, while a translational flex of an infinite frame- 
work is a bounded flex but is not a vanishing flex. Adopting a term 
that has been used in applications [2] we refer to flexes which are not 
bounded as colossal flexes. 

Definition 2.18. A continuous flex p(t) = (p k (t))^ =1 , (t G [0, 1]) of an 
infinite framework (G,p) in R d is said to be 

(i) a bounded flex if for some M > and every k and t, 

\p k (t) - p k {0)\ < M, 

(ii) a colossal flex if it is not bounded, 

(iii) a vanishing flex if p(t) is a bounded flex and if the maximal 
displacement 

\\pk ~Pk (O)!!^ = sup \pk(t) - Pk(0)\ 
te[o,i] 

tends to zero as k — > oo, 
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(iv) a square-summable flex if 

^WPk-PkWWL < OO, 

k=l 

(v) a summable flex if 

oo 

lb* - P*(°)lloo < OO, 

k=l 

(vi) an internal flex if for all but finitely many k the function Pk(t) 
is constant. 

Also we say that (G,p) has a deformation (resp. bounded or vanish- 
ing deformation) if it has a base-fixed flex p(t) (which is bounded or 
vanishing) . 

Definition 2.19. A connected infinite locally finite proper framework 
in two or three dimensions is boundedly rigid (resp. summably rigid, 
square-summably rigid, smoothly rigid, internally rigid) if there is no 
deformation, that is, no base-fixed proper continuous flex, which is 
bounded (resp. summable, square-summable, smooth, internal). 

2.5. Sufficient conditions for flexibility. There is a sense in which 
vanishing flexibility is almost the only obstacle to the existence of a 
flex of a framework all of whose finite subframeworks are flexible. More 
precisely, in the hypotheses of the next theorem we assume that there 
are two distinguished framework vertices, p±,p 2 , such that any finite 
subframework (H,p) containing p±,p 2 has a flex which properly sep- 
arates this pair in the sense of condition (ii) below. The additional 
requirement needed is that there is a family of flexes of the finite sub- 
frameworks whose restrictions to any given subframework (H,p) are 
uniformly smooth in the sense of condition (i). Note that the con- 
stant here depends only on (H, p) and indeed the resulting flex may of 
necessity be a colossal flex. 

Theorem 2.20. Let (G,p) be an infinite locally finite framework in IR d 
with a connected graph, let 

be subframeworks, determined by finite subgraphs G r = (V r , E r ) with 
union equal to G and let t>i,t>2 be vertices in G\. Suppose moreover 
that for each r = 1,2, ... , there is a base-fixed smooth flex p T (t) = 

(P*(*))L=? r)l ofQ r = (G r ,p) such that 

(i) for each finite framework Q\ the set T\ of restriction flexes 



{p r (t)\g r .r>i} 
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have uniformly bounded derivatives, that is, there are constants Mi, I = 
1,2, ... , such that 

\lpl(t)\<M l forr>l,v k eV h 

(ii) the framework points p\,P2 are uniformly separated by each flex 
p r (t) in the sense that 

b r i(i)-p r 2 (i)l-K(o)-^(o)l>c 

for some positive constant c. 
Then (G,p) has a deformation. 

Note that it is essential that the separated vertices of condition (ii) 
are the same for each subgraph. To see this note that the two-way 
infinite trapezium strip framework considered in Figure 7 has smooth 
deformations on each of its finite strip subframe works, each of which 
"separates" some two vertices (at the end of the strip) by a fixed posi- 
tive distance. Nevertheless the infinite strip fails to have a deformation. 

Proof. For / = 1,2, ... , let X\ be the space of continuous functions 
from [0, 1] to and note that the family Ti, by the hypotheses, is 

an equicontinuous family in X;. Moreover with respect to the supre- 
mum norm J 7 ; is a bounded set. By the Ascoli-Arzela theorem (see 
[5U] for example) J 7 ; is precompact and in particular for T\ there is a 
subsequence r 1 ,r 2 , ■■■ such that the restrictions p Tk \Q\ converge uni- 
formly in their (finitely many) coordinates to a flex q 1 of Q\. Relabel 
the sequence (p rk ) as {p^' 1 '). The restrictions of these flexes to Q 2 sim- 
ilarly have a convergent subsequence, say (p^' 2 ^), and so on. From this 
construction select the diagonal subsequence (p( fc ' fc )). This converges 
coordinatewise uniformly to a coordinatewise continuous function 

q : [0, l^l^x^x .... 

Since the restriction of q to every finite subframework is a flex, q satisfies 
the requirements of a flex of (G,p), except possibly the properness 
requirement (v) of Definition 2.1. In view of (ii) however, q(0) ^ q(l) 
and so q is a continuous flex of (G,p). □ 

Remark 2.21. Computer simulations provide evidence for the fact 
that small random perturbations of Q%i yield frameworks that are flex- 
ible. That is it seems that Qj? is stably flexible. It would be interesting 
if the theorem above could assist in a proof of this. 

2.6. Flex extensions and generic rigidity. Let us recall a version 
of Laman's theorem. 

Theorem 2.22. Let {G,p) be an algebraically generic finite framework. 
Then (G,p) is infinitesimally rigid if and only if the graph G has a 
vertex induced subgraph H , with V(H) = V(G), which is maximally 
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independent in the sense that 2\V(H)\ = \E(H) \ + 3 and 2\V(H')\ > 
\E(H') \ + 3 for every subgraph H' of H . 

For convenience we refer to a maximally independent finite graph as 
a Laman graph. We remark that any Laman graph can be obtained 
from a triangle graph by a sequence of moves known as Henneberg 
moves. The first of these adds a new vertex with two connecting edges 
while the second breaks an edge into two at a new vertex which is then 
connected by a new edge to another point of the graph. 

Now let G be an infinite graph which contains a subgraph H on all 
the vertices of G and suppose that H is a-Laman. In view of Laman's 
theorem every algebraically generic realisation of H (and hence G) in 
the plane is u-rigid and so continuously rigid. Is the converse true ? 
That is, if every generic realisation of an infinite graph G is rigid does 
G necessarily contain a u-Laman subgraph with the same vertex set. 

To see this one needs to show is that if G is cr-(Laman-l) and not 
a-Laman, then there exists a vertex generic realisation (G,p) in the 
plane which has a continuous flex. That is, we want to build up a flex 
of the infinite structure by adding new vertices and edges, in the least 
handicapping way, to allow all, or most of the flex of an initial finite 
subgraph to be extended. 

Alternatively, and more explicitly, suppose that one starts with a 
generic connected Laman-1 framework (G\,p) with n vertices and an 
infinite sequence of Henneberg move "instructions". These instructions 
yield a unique infinite graph. Is it possible to choose associated frame- 
work points Pn,Pn+i, ■■ so judiciously that some (perhaps small) flex of 
(Gi,p) extends fully to each successive finite extension framework (and 
hence the infinite framework) ? 

3. Rigidity Operators and Infinitesimal Rigidity 

In previous sections we have considered some variety in the nature of 
continuous flexes p(t) and how they might distinguished. A companion 
consideration is the analysis of various spaces of infinitesimal flexes. 
This gives insight into continuous flexes since the derivative p'(0) of a 
differentiable flex p(t) is an infinitesimal flex. 

Here we give an operator theory perspective for an infinitesimal the- 
ory of infinite frameworks in which the rigidity matrix R(G, p) is viewed 
as a linear transformation or linear operator between various spaces. 
The domain space contains a space of infinitesimal flexes, which lie in 
the kernel of the rigidity operator, while the range space contains a 
space of self-stress vectors namely those in the kernel of the transpose 
of R(G,p). 

3.1. Infinitesimal rigidity and the rigidity matrix. Recall that 
for a finite framework (G,p) in M. d with n = \V\ an infinitesimal flex 
is a vector u = (%,..., u n ) in the vector space % v = M. d © • • • © R d 
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such that the orthogonality relation (pi — Pj,Ui — Uj) = holds for 
each edge (vi,Vj). This condition ensures that if each pj is perturbed 
to Pi(t) = Pi + tUi, with t small, then the edge length perturbations are 
of second order only as t tends to zero. That is, for all edges, 

\Pi(t)-pM\-\Pi-pj\ = o(t 2 ). 

If q(t) : [—1, 1] — > H v is a two-sided smooth flex of the finite frame- 
work (G,p) then q'(0) is an infinitesimal flex and for a generic finite 
framework every infinitesimal flex arises in this way. See Asimow and 
Roth pQ for example. 

Associate with an infinite framework (G,p) the product vector space 

H v = Y[R d = R d ®R d ® ... 
v 

consisting of all sequences u = (ui, u 2 , ■ ■ ■)■ Conceptually such a vec- 
tor corresponds to a specification of instantaneous velocities, or to a 
perturbation sequence, applied to the framework joints. Define an in- 
finitesimal flex of (G,p) to be a vector u in R v for which, as above, 
(pi —pj, Ui — uj) = holds for each edge (vi, Vj), and let Rfi denote the 
linear space of all these vectors. In the planar case Rfi contains the 
three-dimensional linear subspace (assuming G has at least one edge) 
of the infinitesimal flexes that arise from the isometries of M. 2 . Note 
that a nonzero rotation infinitesimal flex u is an unbounded sequence 
if and only if (G, p) is an unbounded framework. We denote the space 
of rigid body motion infinitesimal flexes as R r ig- 

Definition 3.1. An infinite framework (G,p) is infinitesimally rigid if 
every infinitesimal flex is a rigid body motion infinitesimal flex. 

If G is infinite then % v contains properly the direct sum space 
"H°° = Y2v consisting of vectors whose coordinates are finitely 
supported, in the sense of being finitely nonzero. The following defini- 
tion is convenient and evocative. 

Definition 3.2. An infinite framework (G,p) is internally infinitesi- 
mally rigid if every finitely supported infinitesimal flex is the zero flex. 

We now give the usual direct definition of the rigidity matrix R(G,p) 
of a framework (G,p), allowing G to be infinite. This matrix could also 
be introduced via the Jacobian of the equation system that defines 
V(G,p) since 2R(G,p) is the Jacobian evaluated at p. 

Write pi = (xi, yi), Ui = (uf, uf),i = 1,2,..., and denote the coordi- 
nate difference Xi — Xj by Xij . The rigidity matrix is an infinite matrix 
R(G,p) with rows indexed by edges e±, e%, . . . and columns labeled by 

vertices but with multiplicity two, namely vf, vf, v%, v%, Note that 

any matrix of this shape, with finitely many nonzero entries in each 
row, provides a linear transformation from R v to R e = YIe ^- 
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Figure 14. Unit cell for an internally infinitesimally 
flexible periodic framework. 

Definition 3.3. The rigidity matrix of the infinite framework (G,p), 
with p = (pi) = (xi, Hi), is the matrix R(G, p) with entries Xij,Xji, y^, yji 
occurring in the row with label e = (vi,Vj) with the respective column 
labels vf,Vj,v%,Vj, and with zero entries elsewhere. 

In particular a vector u in % v is an infinitesimal flex if and only if 
R(G,p)u = 0. 

Definition 3.4. (i) A stress (or, more properly, a self-stress) of a finite 
or infinite framework (G,p) is a vector w = (w e ) in T-L e = n_B-^ such 
that R(G,pYw = 0. 

(ii) A finite or infinite framework (G,p) is isostatic, or absolutely 
isostatic, if it is innnitesimally rigid and has no nonzero self-stresses. 

Since it is understood here that G is a locally finite graph a rigidity 
matrix hase finitely many entries in each column and so its transpose 
corresponds to a linear transformations from T-t e to r H. v . 

In the finite self-stress represents a finite linear dependence 

between the rows of the rigidity matrix, which one might abbreviate, 
with language abuse, by saying that the corresponding edges of the 
framework are linearly dependent. A self-stress vector w = (w e ) e£ E 
can be simply related to a vector b = (b e ) conceived of as a sequence 
b = (b e ) of bar tension forces with a resolution, or balance, at each 
node. Indeed, for such a force vector b the vector w for which w e = 
\pi — pj\~ x b e (e = (vi,Vj)) is a stress vector. Thus there is a simple 
linear relationship between the space of internal stresses and the space 
of resolving bar tensions. We shall not consider here the more general 
stress vectors, important in engineering applications, that arise from 
an external loading vectors. 

Let "Hg° be the space of finitely supported vectors in "H e . We say 
that an infinite framework (G,p) is finitely isostatic if it is internally 
innnitesimally rigid and if the finite support stress space %°^ r := %str^ 
"Hg is equal to {0}. It is straightforward to see that the grid frameworks 
G%d, in their ambient spaces, are finitely isostatic, as is the kagome 
framework. 
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Between the extremes of infinitesimal rigidity and internal rigidity 
there are other natural forms of rigidity such as those given in the 
following definition. Write £°°, £ 2 and Cq to indicate the usual Banach 
sequence spaces for countable coordinates, and write "H^ ,?^ , . . . ,H® 
for the corresponding subspaces of R e and % v . 

Definition 3.5. An infinite framework (G,p) is 

(i) square- summably infinitesimally rigid (or infinitesimally £ 2 -rigid) 

if 

U 2 V (1 ker R(G,p) = {0}, 

(ii) boundedly infinitesimally rigid (or infinitesimally £°°-rigid) if 

n ker R(G, p) = ^n Rri g , 

(iii) vanishingly infinitesimally rigid (or infinitesially Co-rigid) if 

Ul n ker R(G,p) = {0}, 

(iv) square-summably isostatic (or £ 2 -isostatic) if it is infinitesimally 
£ 2 -rigid and 

n str nn 2 e = {o}, 

(v) boundedly isostatic if it is boundedly infinitesimally rigid and 

UstrMif = {0}, 

(vi) vanishingly isostatic if it is vanishingly infinitesimally rigid and 

n str nn e = {o}. 

There is companion terminology for flexes and stresses. Thus we 
refer to vectors in R? v D ker R(G,p) as square summable infinitesimal 
flexes and so on. 

Example 3.6. Let us use the shorthand (N,p) to denote a semi-infinite 
framework in M 2 whose abstract graph is a tree with a single branch, 
with edges (vi, v 2 ), (v 2 , v 3 ), . . . and where p = (pi),Pi = (xi,yi),i = 
1,2,... . Then, writing Xij and jjij for the differences Xi — Xj and 
Hi = jjj, as before, the rigidity matrix with respect to the natural 
ordered bases takes the form 



R(N,p) 



Xl2 Vl2 X 2 l ?/21 ... 
223 2/23 ^32 2/32 
* * * 



With respect to the coordinate decomposition R v = R x © R y we have 
R(N,p)= [R x R y ] = [D x D y ] 
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where R x 
trices 



Z12 
X 23 

x 34 



DyT, where D x and D y are the diagonal ma- 



D, 



'2/i2 

J/23 

2/34 



and 



T 



I -1 
1 




If we now identify the domain and range spaces in the natural way 
for these coordinates then we have T = I — U t where U l is the trans- 
pose of the forward unilateral shift operator on the linear space of real 
sequences. 

The analogous framework (Z,p) has a similar matrix structure in 
all respects except that in place of the Toeplitz matrix T one has the 
corresponding two-way infinite Laurent matrix / — W~ x where W is 
the forward bilateral shift. In both cases, the two-dimensional sub- 
space spanned by the translation flexes is evident, being spanned by 
the constant vectors in % x and % y . Evidently there are infinitely many 
finitely supported flexes and in fact it is possible to identify ker R(G,p) 
as a direct product space. 

One can use operator formalism to examine the space of stresses. In 
the case of the simple framework (Z,p) note that W l = W~ Y and 

ker(7 -W) = ker(7 - W' 1 ) = Re 

where e is the vector with every entry equal to 1. Since 

\I-W 

I-W 



R(z, P y 





~D X ~ 




D y _ 



it follows that a vector w is a stress vector if and only if D x w G Me 
and D y w G Me. Thus for some constants a, (3 we have x^+iWj = 
a, yi t i + iWi = (3, and so for all i 

y% - Vi+i = P 

This colinearity condition shows that the space of stresses is trivial 
unless the framework points pj, i 6 Z are colinear in which case T-L s tr 
is one dimensional. This includes the colinear cases in which p is a 
bounded sequence and the framework lies in a finite line segment in 



Example 3.7. With similar notational economy write (Z r ,p), (resp. 
(W,p)) for frameworks associated with the grid graph with vertex set 
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labeled by r-tuples of integers (resp. positive integers) n = . . . , n r ) 
where the edges correspond to vertex pairs (n, n±ej), where e±, . . . ,e r 
are the usual basis elements. The ambient space for the framework is 
either understood or revealed by the entries of the vector p. 

Again one can use operator formalism to analyse the space of stresses 
as a vector subspace of H. In the special case of the regular grid frame- 
work Q12 in M 2 one can see that the vector subspace l-i s tr is a direct 
product vector space (like H itself) whose product basis is indexed 
by (two-way infinite) linear subframeworks parallel to the coordinate 
axes. This is also true for a general orthogonal grid framework such 
as the bounded grid framework determined by the framework points 
(±(1- (1/2)*), ±(1 "(1/2)0). 

Note that we have defined an infinitesimal flex in a local way, being 
the verbatim counterpart of the usage for finite frameworks. In par- 
ticular the notion takes no account of the possibility of (second order) 
amplification or vanishing flexibility. 

Our examples above indicate the importance of shift operators and 
in the next section we see that the bilateral shift operators, in their 
Fourier transform realisation as multiplication operators, play a central 
role in the discussion of periodic frameworks. 

Let us note here that approximate infinitesimal flexes are natural for 
infinite frameworks and the operator theoretic perspective. 

Definition 3.8. An approximate square-summable flex of an infinite 
framework (G, p) is a sequence of finitely supported unit vectors U\, u 2 , ■ ■ 
in "H 2 (or /C 2 ) such that ||i?(G,p)M n || 2 — > as n — y oo. 

Let (G,p) be a distance regular framework. Then it is straightfor- 
ward to show that the rigidity matrix determines a bounded Hilbert 
space operator R. It is the metrical and geometric properties of the ac- 
tion of R and its transpose that have relevance to rigidity theory rather 
than the spectral theory of R. However we do have the almost vac- 
uous statement that the existence of approximate (square-summable) 
flexes corresponds to the point belonging to the approximate point 
spectrum of R. For if lies in the approximate point spectrum then 
(by definition) Rv n is a null sequence for some sequence of unit vec- 
tors v n , and approximation of these unit vectors by vectors with finite 
support yields, after normalisations, an approximate square-summable 
flex sequence («„). 

It is implicit in the matricial function association below that the 
rigid unit modes of translationally periodic frameworks are tied to the 
existence of approximate flexes. 

4. Crystal Frameworks and Flexibility 

In previous sections we have constructed frameworks to illustrate 
various definitions and properties. It is perhaps of wider interest to 
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understand, on the other hand, how extant infinite frameworks, such 
as those suggested by crystals or repetitive structures, may be flexible. 
Accordingly we now define crystal frameworks and investigate various 
forms of flexibility and rigidity 

4.1. Periodic and crystal frameworks. We have already observed 
some properties of the basic examples of the grid framework Qj2 , the 
squares framework Q sq and the kagome framework, Qkag- In we 
also have analogues, such as the cube framework Q cu be, the octagon 
framework Q oct and the kagome net framework Qknet, which consist, 
respectively, of vertex-joined cubes, octahedra and tetrahedra, with no 
shared edges or faces, each in a natural periodic arrangement. The 
frameworks Qkag, Qcube, Qoct and Qknet are polytope body-pin frame- 
works but we consider the polytope rigid units as bar-joint subframe- 
works formed by adding some, or perhaps all, internal edges. As such 
these frameworks are examples of crystal frameworks in the sense of 
the formal definition below. We first comment on a wider notion of 
periodicity. 

Definition 4.1. An affinely periodic framework in M d is a framework 
Q = (G,p) for which there exists a non-trivial discrete group of affine 
transformations T g ,g G V, where each T g acts on framework points 
and framework edges. 

For example, the two-way infinite dyadic cobweb framework of Figure 
15 is affinely periodic for the dilation doubling map and the four-fold 
dihedral group D4. 




Figure 15. An affinely periodic cobweb framework. 



For another example we may take the infinite Z-periodic framework 
in three dimensions for which Figure 16 forms a perspective view down 
a central axis, with framework vertices (±1,±1, m),m G Z. Here the 
affine group is an isometry group isomorphic, as a group, to Z x C2 x D^. 
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To illustrate the following definition observe in Figure 16 a template 
of six edges and three vertices which generates the kagome framework 
by the translations associated with the parallelogram unit cell. Bor- 
rowing crystallographic terminology we refer to such a template as a 
motif for the framework and the chosen translation group. The fol- 
lowing formal definition gives a convenient way of specifying abstract 
crystal frameworks. 



Figure 16. A motif and unit cell for the kagome framework. 



Definition 4.2. A crystal framework C = (G,p) in M. d is a connected 
bar-joint framework for which there is a discrete group of translation 
isometries T = {T g : g £ T>}, a finite connected set F e of frameworks 
edges and a finite set F v of framework vertices (being a subset of the 
vertex set of F e ) such that 
(i) the unions 



coincide with the sets of framework edges and vertices, and 
(ii) these unions are unions of disjoint sets. 

We denote a crystal framework C by the triple (F v , F e , T) or by the 
triple (F vl F ej Z d ) in the case of integer translations. An associated 
unit cell for C may be defined as a set which contains F v and no other 
framework points and for which the translates under the translation 
isometries are disjoint and partition the ambient space. For example 
in the case of Qy2 we may take the semiopen set [0, l) 2 or the set 
[0, 1) x [1/2, 5/2) as unit cells. Such parallelepiped unit cells are useful 
for us for torus models for crystal frameworks. Voronoi cells (Brillouin 
zones) also play a unit cell role in applications but we shall not need 
such geometric detail here. 

In many applied settings the appropriate framework models have 
"short" edges, spanning no more than two adjacent unit cells. Here we 
allow general edges which may span a chain of adjacent cells. 

Recall from elementary crystallography that, modulo orthogonal trans- 
formations, there are 14 different forms (or symmetry types) in which 
a countable set of isolated points can be arranged with translational 
symmetry throughout three-dimensional space. These arrangements 




^gevT g (F e ) 



UgeT>T g (F v ) 
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are called the Bravais lattices and the translation group T above cor- 
responds to such a lattice. Thus each point of the framework lies in 
the Bravais lattice generated by the orbit of its unique corresponding 
motif vertex under the translational group. 

4.2. Deformability and flow flexibility. Recall that a general (count- 
able, locally finite, connected) framework Q is rigid if there is no 
base-fixed continuous flex and is boundedly rigid, or boundedly non- 
deformable, if there is no bounded base- fixed continuous flex pit)- Re- 
call, from Section 2, that bounded flexes are those for which there is an 
absolute constant M such that for every vertex v the time separation 
\Pv{t) — Pv(0)\ is bounded by M for all t and all v. 

We first describe a context for the standard "alternation" flexes of 
Q aq and Q kag and certain periodic flexes of Qkag with reduced symmetry. 
For these nonbounded flexes translational periodicity is maintained but 
relative to an affine flow of the ambient space. By an affine flow we 
mean simply a path t — > A t of affine transformations of M. d which is 
pointwise continuous. The simplest such flow in two dimensions is a 
contracting flow such as A t (x, y) = ((1 — tc)x, (1 — tc))y, < c < 1. 
The alternating flexes of Q sq and Qkag are associated with such a flow. 

Definition 4.3. Let d = 2,3, let t A t be a flow of R d and let C = 

(G, p) be a crystal framework for the translation group T = {T g : g G 
V}. A flow-periodic flex of C, relative to the flow and the translation 
group T, is a continuous flex p(t) such that for each t the framework 
(G,p(t)) is V- periodic, where T = {A t T g A~l x : g G T>}. 

Let us say simply that C is affinely periodically deformable if there 
exists a non trivial flow-periodic flex. 

A flow-periodic flex p{t) for C can be defined for a given flow if and 
only if for each t one can continuously solve the distance constraint 
equations for the vertex positions of the motif, with the periodicity 
constraint, in the A t - deformed unit cell. Equivalently, the motif and 
unit cell define a finite framework on a torus with the noninterior edges 
of the motif providing reentrant ("locally geodesic") edges on the torus. 
For example consider the torus framework in Figure 17. 




FIGURE 17. A torus framework. 
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A horizontal affine contraction H t : (x,y) —> ((1 — t)x,y) leads to a 
continuous flex of the torus framework and hence to a colossal flex of 
the associate crystal framework. The same is true for the vertical affine 
contraction Vj and for the skew affine transformation 

S t : (x, y) — > (x + (smt)y, y + (1 — cost)x). 

Note that this particular transformation preserves both the cyclic width 
and cyclic height of the torus. 

For an illustration of the method we note the following affine de- 
formation result for what might be termed periodic cell-generic grid 
frameworks in ~R d . For more general results of this nature see also 
Borcea and Streinu [2]. 

Theorem 4.4. Let C = (1* d ,p) be a grid framework in M. d which is 
(rii, . . . , rid) -periodic o,nd which is a 1/3-perturbation ofQ^d i n the sense 
that 

\p(kx, . . . , k d ) - . . . , k d )\ < 1/3 

for < ki < ni — 1,1 < i < d. Then C is affinely periodically de- 
formable. 

Proof. We sketch the proof in case d = 2 with (ni,ri2) = (n,m). For 
convenience re-scale the framework so that the large cell [0, n — 1) x 
[0,m — 1) of C becomes to the usual unit cell. Let (F v ,F e ,Z 2 ) be the 
motif representation of (Z 2 ,p). Suppose first that n = m = 2 and the 
motif has four noninterior edges (corresponding to the four reentrant 
edges on the associated torus). Let F' e be F e with one vertical and 
one horizontal reentrant edge removed and let C be the associated 
framework. (See Figure 17.) For t taking positive values in some finite 
interval there are affine deformations of C associated with each of the 
flows H t ,Vt,St- Moreover there is an affine deformation associated with 
any composition A t = Vp(t)H a (t)Su where the functions a and (3 are 
any continuous functions with a(0) = /3(0) = 0, where t takes values 
in some finite interval. Now note that we may chose a(t) so that the 
separation distance corresponding to the omitted horizontal edge is 
constant, and we may choose /3(t) similarly for the omitted vertical 
edge. Note also the essential fact that VpM does not change the cyclic 
width. In this way an affine deformation for C is determined. (In the 
case of generic points, is determined uniquely on some finite interval). 

The same principle operates for general (n,m). The subframework 
C, with two deleted noninterior edges (one for each coordinate direc- 
tion) once again has a flow periodic deformation associated with a com- 
position At = V/3(t)H a (t)St, for t in a sufficiently small interval. This 
follows from a simple induction argument. Once again the functions 
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a and /3 are chosen to provide a deformation of C in which the sepa- 
rations for the omitted edges is constant. In this way a flow periodic 
deformation of C itself is determined. □ 

The deformations obtained in the last proposition are colossal de- 
formations and this seems to be a necessary condition if some form of 
periodicity is to be maintained. It would be interesting to determine 
when such frameworks possess bounded "unstructured" deformations. 

One may also consider affine deformations relative to a subgroup of 
T> associated with a supercell. By a supercell we mean a finite union 
(or possibly an infinite union) of adjacent cells which tile the ambient 
space by translations from the subgroup. In this case the deformations 
maintain only a longer period of translational symmetry and for an 
infinite linear supercell one may even forgo translational symmetry in 
one direction. 

To illustrate this let us note a class of interesting deformations of the 
kagome framework (G,p) which have this form. 

Start with an alternation flex p(t) which leaves fixed a particular 
vertex and leaves invariant a line of hexagon diameters. Note that 
this flex (which is not base-fixed) has bilateral symmetry with respect 
to the line (viewed as a mirror line) but that the inversion symmetry 
of {G,p) about the fixed vertex is broken for t > 0. Perform identical 
surgeries on the frameworks (G,p(t)) as follows: 

(i) cut the frameworks along the fixed line, 

(ii) effect a reflection of one of the resulting half-planes frameworks 
in the orthogonal line through p* and rejoin the half-planes to create a 
new flex of (G,p(0)). 

We might view this flex as one with a symmetry transition line. The 
hexagons in this line maintain inversion symmetry while elsewhere the 
hexagons maintain bilateral symmetry. 

One may perform such surgery on several parallel surgery lines si- 
multaneously. Performing surgery on countably many lines leads to 
the following theorem. This shows, roughly speaking, that the kagome 
framework sits in its configuration space as an infinitely singular point 
in the sense that it is the starting point for uncountably many distinct 
flexes. (Here we use the term distinct in the sense of Definition 2.7.) 

Theorem 4.5. There are uncountably many distinct flow periodic flexes 
of the kagome framework. 

4.3. Crystal frameworks and periodic infinitesimal flexibility. 

Let C = (G,p) = (F v , F e ,7* 2 ) be a crystal framework in the plane for 
the integer translation isometry group. A natural form of infinitesimal 
flex for C is that of a 1-cell periodic flex in the following sense. 

Definition 4.6. A vector u = (u v ) v& v in the real vector space H v , or 
in the complex vector space K. v , is a 1-cell-periodic infinitesimal flex for 
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C (or simply a periodic flex if there is no ambiguity) if u G ker R(G,p) 
and u K+n = u K for all k G F v and n G Z d . 

Such an infinitesimal flex a is a bounded sequence of vectors deter- 
mined by periodic extension of what we may call the motif flex vector 
Umotif = (u K ) Ke F v . There is a one-to-one correspondence between these 
periodic flexes and the finite vectors that are in the kernel of the motif 
rigidity matrix R m (G,p) which we may define as the natural "periodic 
completion" of the \F e \ xd\F v \ submatrix of R(G,p). This is the natural 
representing matrix for R(G,p) viewed as a linear transformation be- 
tween the finite dimensional subspaces of Z 2 -periodic vectors. If there 
are no reflexive edges in the motif then 

Similarly we define periodic infinitesimal stresses as those which cor- 
respond to periodic extensions of vectors in the cokernel of the motif 
rigidity matrix. 

The crystal framework C, with given motif and discrete translation 
group indexed by V, is said to be 1- cell-periodically isostatic (or simply 
periodically isostatic) if the only 1-cell periodic flexes are translation 
flexes and if there are no nontrivial 1-cell periodic stresses. 

Once again it can be helpful to consider a flat torus model for such 
frameworks. Note for example that the generic periodic framework de- 
fined by Figure 17 is periodically rigid and indeed periodically isostatic. 

The following interesting periodic variant of Laman's theorem has 
been obtained recently by E. Ross [28]. Let us say that the planar 
periodic framework C is topologically proper if C is connected and for 
every torus sub framework motif F^, F' e with 2|i^| — |.Fg| =2 there is an 
edge cycle from F' e which properly wraps around the torus in the sense 
that the associated homotopy class is nonzero. 

Theorem 4.7. Let C = (F v , _F e ,Z 2 ) be an infinite framework in the 
plane which is periodic for the integer translation group and is topolog- 
ically proper. Then following are equivalent. 

(i) 2\F V \ — \F e \ = 2 and for all edge induced submotifs F^, F' e we have 
2\F' v \-\F' e \>2. 

(ii) C is periodically isostatic in the sense that the periodic vectors in 
ker R(G,p) are spanned by the two translations (periodic rigidity) and 
the periodic vectors in coker R(G,p) are zero (periodic stress free). 

As we have seen, a periodically infinitesimally rigid framework in 
the sense above may nevertheless correspond to a flat torus framework 
which can flex infinitesimally (and even deform) relative to contraction 
or expansion of the unit cell. Allowing such vertical and horizontal 
freedoms adds two further degrees of freedom to the constraint equa- 
tions. It would be of interest to obtain a similar characterisation in 
this well as for the triply flexible affine torus. 
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5. The Matricial Symbol Function of a Crystal 

Framework 

We now derive matrix function operators for general (abstract) crys- 
tal frameworks. Our approach is decidedly Hilbert space theoretic 
and allows for an extended conceptual framework for rigidity analy- 
sis. First we outline the standard identification of operators commuting 
with shift operators and multiplication operators in a Fourier transform 
space. We adopt complex scalars, replacing the % v spaces of Section 3 
by their complex counterparts, K v ,K, e , K, 2 , and so forth. 

5.1. Matrix function multiplication operators. Let be a con- 
tinuous complex- valued function on the unit circle T. It defines a mul- 
tiplication operator T on the usual complex Hilbert space L 2 (T). Its 
representing matrix with respect to the standard orthonormal basis 
{z n : n G Z} is the ZxZ indexed matrix with entries 

T hJ = (^,z l ) = (^-\l) = 4>( i -j), 

the (i — j) th Fourier coefficient of 4>. 

Similarly, let <&(z) be a continuous matrix- valued function on the 
two or three dimensional torus T d taking values in the space of n x m 
complex matrices. One can specify such a function = <3>(zi, . . . , zj) 
in terms of a matrix of scalar functions $(z) = [4>k,i{ z )]k'^i i=i- Given 
in this way <&(z) defines a multiplication operator from L 2 (T) cg> C m to 
L 2 (T) ®C n . Indeed, choose a basis {^} for C m and an associated basis 
{6 <S> z p } for the domain space. Similarly let {r]k <g> z p } be a basis for 
the codomain. Then the operator T of multiplication by &(z) may be 
defined by specifying, for each continuous function f(z) on T d , 

in 

k=i 

and extending by linearity and continuity. Thus we can see that the 
representing matrix for T with respect to these bases is determined by 
the Fourier coefficients of the matrix entries for $: 

(T(& ® z p ),T] k <g> z q ) cm ® L 2 {Td) = ((f) k j(z)z p - q , i) L2(Td) = 4>k,i(q-p)- 

Viewing an element of L 2 (T d ) ® C m as a function F(z) taking values in 
C m (strictly speaking, taking values almost everywhere), and similarly 
for vectors in the codomain, one also considers TF as the function 
$(z)F(z). This operator is usually denoted as M$. 

Suppose now, starting afresh, we have orthonormal bases 

: 1 < I < m,p G Z d } { VktP :l<k<n,pe Z d } 

for Hilbert spaces such as K 2 and K 2 respectively. Suppose also that 
we are given an operator T from K 2 to K 2 by means of its representing 
matrix with respect to these bases and suppose moreover that this 
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matrix has the translational symmetry above in the sense that the 
matrix entries 

(T£i, p , Vk,q)tc 2 e 

are independent of p — q. Let F v : C n (g) L 2 (Y d ) — > /C 2 be the Fourier 
transform and T e the Fourier transform for /C 2 . These are simply the 
unitary operators determined by the natural bijection of basis elements, 
namely — > £1 <S> z p and r}k, P —> f]k ® z p , respectively. Thus, from 
the scalar matrix entries the functions <pk.i are determined and hence 
the matrix function Q(z). In this way one identifies F~ X TF V as the 
associated multiplication operator _M$. 

5.2. The rigidity matrix as a multiplication operator. Let us 

keep in view the special case of the Z 2 -periodic grid framework deter- 
mined by an algebraically generic quadrilateral in the unit cell. 

Let p l , . . . ,p A be four framework vertices in the unit cell [0, l) 2 con- 
stituting the motif set F v of Z 2 -periodic quadrilateral grid (Z 2 ,p). Let 
Pij — V t + {h j) j f° r ihj) £ be the general framework points and let 
e 1 , . . . , e 8 denote the eight edges which form the motif set F e given by 

e 1 = \p\p\e 2 = [p\p%e 3 = \p 3 ,p%e 4 = \p\p\ 

e 5 = [p*i,o,PV 6 = K_!,p 2 ],e 7 = K ,p 3 ],e 8 = [p l ^p% 

and let e\ • = e* + be a typical edge. Write the corresponding 
basis for the vector space KP e as 

{te:l<t<8,(z,j)e^ 2 }- 

Also, label the natural basis for JC® as the set of vectors ., Cs,i,j}j 
where s ranges from 1 to 4 and range through Z 2 . 

The rigidity matrix R(G,p) satisfies the following symmetry equa- 
tions, 

W iJ R(G,p)U itj =R(G,p) 

with respect to the translation shift operators Uij on /C°° and Wij on 
JC®°. By the discussion above the rigidity matrix for this framework 
gives rise to an 8 x 8 matrix function $(^1,22) = (<f>i j(zi, z 2 ) on the 
torus T 2 in C 2 which is determined by the equations 

together with a companion set of equations for the ^/-labeled basis el- 
ements. Furthermore there is a simple algorithm for computing the 
symbol matrix function from the motif (F v , F e ) and the motif rigidity 
matrix. 

For a general crystal frameworks in the plane we have the following 
recipe for identifying $(^1,^2). There is an entirely similar identifica- 
tion of the matrix function of a crystal framework in higher dimensions. 
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Theorem 5.1. Let be the matricial symbol function of the crystal 
framework C = (F v , _F e ,Z 2 ). The entry of the rigidity matrix deter- 
mined by the edge of F e and the column labeled • (resp. v y si ^) 
provides the {—i,—j)-th Fourier coefficient of 4>k,{x,s] (resp. <fik,{y,s})- 

Note that for the motif for the periodic quadrilateral grid each edge 
has vertices with different s index. This is rather typical and in such 
cases it follows that <pk,{x,i} is either zero or has one nonzero Fourier 
coefficient. In the case of a "reflexive" edge in F e , of the form \p K> o,p Kt s\ 
there are entries in the columns for k corresponding to (l-z s )v e where 
v e is the usual vector of coordinate differences appearing in the rigidity 
matrix for the edge e. That is, v e = p Kj o — p K ,s- 

Thus we may obtain the following identification where the direct 
construction of $ from C is given above. 

Theorem 5.2. Let C = (F v , F e , Z d ) be a crystal framework with m = 
\F e \, n = \F V \ and with rigidity operator R from K? v to K? e . 

(i) The motif (F v ,F e ) determines a matrix-valued trigonometric func- 
tion $ : T d — > M mrfn (C) for which there is a unitary equivalence 
F~ 1 RF V = Mq, where M$ is the multiplication operator 

M$ : C nd (8) L 2 (T d ) -> C m (8) L 2 (J d ). 

(ii) The framework C is square- summably rigid (resp. square- summably 
stress-free) if the column rank (resp. row rank) of <&(z) is maximal for 
almost every z G T d . 

The following corollary shows that the existence of a square summa- 
ble flex is a rather strong condition for a crystal framework and in fact 
"typical" frameworks in Maxwell counting equilibrium have no such 
flexes. What is of particular significance for such frameworks however 
is the presence of approximate flexes, in the sense noted in Section 3, as 
these are related to the local rank degeneracies of <&(z) and to periodic 
flexes. 

Corollary 5.3. The following are equivalent for a crystal framework 
with Maxwell counting equilibrium. 

(i) C has an nonzero internal (finitely supported) infinitesimal flex. 

(ii) C has an nonzero summable infinitesimal flex. 

(Hi) C has an nonzero square- summable infinitesimal flex. 

Proof. Note that (i) implies (ii) and (ii) implies (iii) so it remains to 
show (iii) implies (i). 

Suppose that the matrix has a square-summable vector f(z) 
in its kernel. Then det($(z)) = for almost all points in the support 
of f(z). Since det<&(z) is a mutivariable trigonometric polynomial it 
necessarily vanished identically, that is the polynomial det$(z) is the 
zero polynomial. 
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Recall that for any square matrix X over a ring we have 
XX = det(A)/ n where X is given by the usual formula 
Xij = (— l) l+ ^Xij det(JT y ) involving the cofactors In particular 

if X is not zero and det(X) = then for a nonzero column vector f(z) 
of X we have Xf(z) = 0. Applying this to X = $ provides the desired 
polynomial (finitely supported) vector in this case. 

If X happens to be the zero matrix we make use of the minimal 
polynomial lemma below. To apply the lemma let q be the minimum 
polynomial of X so that q(X) = 0. By the lemma we have q(X) = 
Xqi(X) and qi{X) ^ by minimality. Thus there is a nonzero vector 
f(z) in the range of qi(X) and this provides the desired local flex. □ 

Lemma 5.4. Let q(X) in R[X] be the minimal polynomial of an n by n 
matrix X whose entries lie in an integral domain R and suppose that 
detX = 0. Then q(0) = 0. 

Proof. Let p(X) be the characteristic polynomial det(A/ n — X). Recall 
from the Cayley-Hamilton theorem that p(X) = 0. In the algebraic 
closure of the field of fractions of R the linear factors of q(X) and p(X) 
agree. See Lange [22] for example. Thus 

q(X) = l[(X-a i r, p(A) = H(X - 

with Si < Ti for each i. In particular for suitably large n the polynomial 
q(0) n is divisible by p(0). Since detX = it follows that p(0) = and 
so q(0) n = 0. Since q(0) is in R it is equal to zero and the proof is 
complete. □ 

5.3. Prom motifs to matrix functions. We now consider some ex- 
amples. Adjusting notation, write z, w for the coordinates variables of 
T 2 so that a trigonometric polynomial takes the form a finite sum 

(j)(z,w) = ^0(2,J>*W J . 

hi 

Examining the recipe of Theorem 5.1 above in the case of the periodic 
quadrilateral grid we see that the internal edges e 1 , . . . , e 4 of the motif 
provide four rows for $(2, w) whose entries are constant functions with 
constants corresponding to those of the rows for the rigidity matrix 
R(G,p). The other edges provide rows according to the following sim- 
plified monomial rule. As we have intimated above, this rule applies 
generally when all the " external" edges of the motif have vertices with 
distinct s index: 

the entries of the e th row of§(z,w) that appear in the columns for 
the external (non-motif) vertex are the corresponding entries for the 
motif rigidity matrix multiplied by the monomial 1?W , where {i,j) is 
the shift index for the cell occupied by the external vertex. 



36 



JOHN OWEN AND STEPHEN POWER 



Thus we obtain the matrix identification in the next proposition. We 
write Pi = (xi,yi), for i = 1, . . . , 4. Note, for example, that the entry 
X12 + 1 = X\ — (x2 — 1) denotes the constant function corresponding to 
the entry of R(G,p) for row e 5 and column x\, while — (a; 12 + l)z is the 
entry corresponding to the x-coordinate of the external vertex of e 5 . 



Proposition 5.5. The matricial symbol function $(2;, w) of the peri- 
odic quadrilateral grid framework (W?,p) for the four vertex motif is the 
matrix function on T 2 given by 
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-Z12 


-2/12 




















^23 


2/23 


-£23 


-2/23 
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-£34 
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X12 + 1 


2/12 
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2/23 + 1 
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-(2/23 + l)w 




















X34 - 1 


y34 


(1 - x 34 )z 


-2/34 Z 


—xnvJ 


(1 - 2/41 )uJ 














X41 


2/41 - 1 



Furthermore, if the four points in the unit cell have algebraically in- 
dependent coordinates then (Z 2 ,p) is square- summably isostatic. 

Proof. The first four rows are linearly independent since the points 
Pi, . . . ,p4 have algebraically independent coordinates. It can be shown 
from elementary linear algebra for almost every z, w the full matrix has 
rank 8. Thus Theorem 5.2 applies. □ 

As a special case we obtain the matrix symbol function for Q1 2 de- 
termined by the eight-edged motif, namely 



§(z,w) = 1/2 
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-1 



The determinant in this case is 

~l^zw(l-zf{l-wf. 

On the other hand for this symmetric quadrilateral framework Qj? 
the formalism above can be applied to a smaller unit cell with the 
simple two-edged motif of Figure 18. 
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Figure 18. A motif and unit cell for Q Z 2. 



Note that the index s for the recipe takes the value 1 only and both 
motif edges have vertices with the same s index. It follows that the 
rigidity operator of Qj? associated with the motif is unitarily equiva- 
lent to the multiplication operator on the Hilbert space C 2 ® L 2 (T 2 ) 
determined by the matricial symbol function 

- 1 



^(z,w) 







w 



One can similarly verify the following. 

Proposition 5.6. (i) The matricial symbol function of the kagome 
framework determined by the six-edge motif is given by 
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V3 


zw 


— y/3zw 
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y/3w 








-1 


-V3 . 



(ii) The determinant of $ka g (z,w) is a multiple of 

zw(z — l)(w — l)(z — w) 
(Hi) The kagome framework is square- summably isostatic. 

We remark that the function matrix association above has also arisen 
in engineering in the analysis of Hutchinson and Fleck [19] of the 
stresses and rigidity of the kagome repetitive truss framework. This 
is derived from the crystallographic perspective of wave periodic flexes 
and Bloch's theorem. 

5.4. Symmetry equations for infinite frameworks. Let C be a 

crystal framework in M d with complex Hilbert space rigidity operator 
R. Also let K-l, /C 2 be as before, let K, 2 ^ be the space of square-summable 
infinitesimal flexes, let JC 2 tr be the space of square-summable infinites- 

K? v e 1C 2 fl , M = /C 2 K? str be their comple- 



imal stresses and let Ai 
mentary spaces. 
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As well as commuting with the coordinate shift operators the rigidity 
operator satisfies commutation relations for every isometric symmetry 
of C. For example, let S be an isometric (not necessarily linear) op- 
erator on M. d which effects a symmetry of C. There is an induced 
unitary operator S e on the complex Hilbert space £ 2 (E) and an analo- 
gous operator S v on £ 2 (V) which permute the standard basis elements. 
Write S v for the isometric (not necessarily linear) operator S v £g> S on 
1~Lv — £ 2 {V) £g> where S is also viewed as an isometric operator on 
C d . Then we have the fundamental symmetry equation S e R = RS V . 
(See JZTJ) 

Moreover, if Q is the symmetry group of C arising from isometric 
transformations of M. d and if p e and p v are the associated representa- 
tions of Q on /Cg and K 2 then we have the symmetry equations 

p e (g)R = Rpv(g) for g e Q. 

We have shown in [27] how such symmetry equations may be used to 
obtain a simple proof of a unitary equivalence which implies the Fowler- 
Guest formula [9] together with various generalisations. In the present 
setting we have the following analogue which also leads to counting 
conditions for isostatic and rigid frameworks. 

Theorem 5.7. Let C be a crystal framework with isometry symmetry 
group Q and let 

Pv = Pm® Pfl 

Pe = PN® Pstr 

be the decompositions associated with the spaces of square- summable 
infinitesimal flexes and stresses. Then pj^ and are unitarily equiv- 
alent representations. In particular if C is square summably isostatic 
then p v and p e are unitarily equivalent. 

Proof. In the square-summably isostatic case R is a bounded operator 
with trivial kernel and trivial cokernel. The partially isometric part U 
of the polar decomposition R = U(R* i?) 1//2 is therefore unitary and it 
is a standard verification that this unitary also intertwines the repre- 
sentations. In general the symmetry equations show that the space JC 2 ^ 
is reducing for p v and K? str is reducing for p e and so the asserted direct 
sum decompositions do exist. Now the restriction of R to M. maps 
to M with trivial kernel and cokernel and so as before the restriction 
representations are unitarily equivalent. □ 

5.5. From matrix function to wave modes. For a crystal frame- 
work C = (F v ,F e ,Z d ) with a given motif let $ : T d ->■ M mM be the 
associated matricial symbol function. 

Definition 5.8. The mode multiplicity function of C associated with 
the given motif and translation group is the function p : T d — > Z + 
given by p(z) = dimker$(2). 
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For d = 2 we also consider /i as being parametrised by coordinates 
s,t in [0,1) x [0,1) so that (z, w) G T 2 corresponds to (e 2ms , e 2mt ). 
From the determinant calculations above we obtain the following. 

Proposition 5.9. (%) For the grid framework and the 8-edged motif 
the mode multiplicity function has values yu(0, 0) = 2, 

0(s,o) = MM) = i, 

if s and t are nonzero, and is zero otherwise. 

(ii) For the kagome framework Qkag and the 6- edged motif the mode 
multiplicity function has values /i(0, 0) = 2, at the origin while [i(s, s) = 
fi(s, 0) = /i(0, t) = 1 if s ^ and t ^ and is zero otherwise. 

Consider now the wave flexes of C which we define as the infinitesimal 
flexes which are 1-cell-periodic modulo a phase factor. 

Definition 5.10. Let C = (F v , E v , Z d ) be a crystal framework in M. d . 
An (infinitesimal) wave flex of C is a complex infinitesimal flex u = 
(u v ) ve v which is wave periodic (or more precisely 1-cell wave periodic) 
in the sense that for some vector q in R d 

for each vertex k in the motif set F v and each d-tup\e n = (ni, . . . , n^). 

The values of the mode multiplicity function corresponds to the di- 
mension of the spaces of wave flexes. To see this suppose that <3>(z) 
is the matricial symbol function for C = (G,p) = (F v , F e , Z 3 ), that 
w G T 3 and that det$(w) = 0. Then Q(w)u m = for some nonzero 
complex motif vector u m . For the Dirac delta function S w (z) on T 3 we 
have, informally, <&(z)F(z) = for all z on the 3-torus where F(z) is 
the function S w (z)u m . Thus, taking Fourier transforms it follows that 
the wave periodic vector u = F(F) in K v satisfies R(G,p)u = 0. This 
shows that the bounded (and non square-summable) vector 

U = (^ K ,n) K eF„,neZ3 = F{F){K,n) = w n u K 

is a wave flex. The Dirac delta argument can be rigourised in the usual 
manner and so we obtain the following theorem. 

Theorem 5.11. Let C = (F v , F v ,Z d ) be a crystal framework in M. d 
with associated symbol function $(z) and suppose that C has no inter- 
nal (finitely supported) self-stresses and no internal (finitely supported) 
flexes. Then 

(i) d\F v \ = \F e \, the associated matricial symbol function <&(z) on 
the d-torus is square with full rank almost everywhere and C is square- 
summably isostatic. 

(ii) infinitesimal wave flexes (that are 1-cell periodic-modulo-phase) 
exist, with phase factor w G T d , if and only if det <&(«;) = 0. In 
this case the dimension of the corresponding space of wave flexes is 
dimker $(iy). 
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It follows also from this that the zero set of the determinant of the 
matricial symbol function associated with the motif of an appropriate 
supercell coincides with the phases of supercell-periodic wave flexes. 

Of particular computational and theoretical interest is what one 
might refer to as the wave flex acquisition when a crystal framework de- 
forms under a colossal flex to a framework with higher symmetry. This 
phenomenon serves as a model for the appearance of so-called Rigid 
Unit Modes (RUMs) in higher symmetry phases of various material 
crystals. 

The following theorem generalises an interesting result of Wegner [32J 
for tetrahedral crystals. It may be viewed as an expression of the simple 
principle that additional symmetry often entails additional flexibility. 
Our proof applies to arbitrary crystal frameworks and is quite direct, 
benefitting somewhat from the economy of operator theory formalism. 
We remark that there are also natural operator algebra perspectives 
that are relevant to symmetry considerations. The result shows in 
particular for d = 3 that the RUM set is typically a union of surfaces, 
being the zero set of a single real-valued polynomial, rather than the 
intersection of the zero sets of the real and imaginary part of a complex 
polynomial. 

For a given crystal framework and motif write Q for the subset of 
T d formed by the phases to of the 1-cell periodic wave flexes. We also 
refer to this informally as the RUM set or the RUMs of C. 

Theorem 5.12. LetC = (F v ,F e ,1, d ) be a crystal framework inM. d with 
d\F v \ = \F e \ and suppose that C possesses inversion symmetry. Then 
the set Q has the form 

n = T d n V(p) 

where V(p) is the zero set of a complex polynomial p(z\, . . . , Zd,z\ . . . ,Zd) 
which is real-valued on T d . 

Proof. Let us denote the set of framework vertices and edges as 
V = {k + n : k e F v ,n e Z d } 

£ = { e + n:eeF v ,neZ d } 

Effecting a translation, if necessary, we may assume that the inversion is 
a : x — > —x, that the unit cell has inversion symmetry and o~(F v ) = F v . 
It may or may not be possible to re-choose the noninternal edges of 
F e so that that c(F e ) = F e . Suppose first that this is the case. We 
show that det 3>(z) itself is either real- valued or purely imaginary, from 
which the stated form for Q follows. 

With the notation of Section 5.4 we have the symmetry equation 

p e (cr)R = Rp v {a). 

Recall that p e (cr) and p v {o~) are the isometric operators induced by 
the (isometric) symmetry element a. Taking Fourier transforms this 
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equation takes the form 

U(a)M Hz) = M Hz) V{a) 

where U (a) and V(a) are the unitary operators determined by their ac- 
tion on the distinguished orthonormal bases. In view of the assumption 
we have 

U(a)(rj k 0z l ) =v* e (k)®z~\ k = 1, . . . , |F e | 

(for certain induced permutations cr e ,a v ) and 

V(a)(C k ' ® z l ) = Q {k) (8) z~\ k = 1, . . . , \F v \,i = 1, . . . , d. 

That is, these operators have the form 

U(a) = E a ®J, V(a) — V a <S> J 

where J is the inversion unitary operator on L 2 {T d ) given by (Jf)(z) = 
f(z), and where E a and V a are scalar permutation matrices induced 
by a. Substituting these forms we see that 

$(*) = (E a ® jy 1 ^)^ ® j) = ^^(^k. 

It follows that 

det = det $(z) det( J B (T )~ 1 det(K) = (-l) r det($(z)) 

for some integer r, since the determinant polynomial has real coeffi- 
cients and so the determinant is either real or purely imaginary. 

In the general case because of edges lying on inversion axes (lines 
through the origin) a cannot be made to act freely on the motif edges. 
However for each motif edge index k there is a monomial shift factor 
z p(k) such that with D the diagonal matrix function 

D = diag{z p[l) , . . . , z pW ) 

we have 

U(a) = D(E a <g) J). 

Now we see that det $(2) = (— l) T 2; :p det( < l ) (z)) for some multi-index p. 
Write $ = F x + iF 2 with the F t real-valued and (-l) r z p = e i7(2:) with 
7 = j(z) real. Equating real and imaginary parts it follows that either 
C0S7 = 1 in which case det<J?(z) = det($(z)) and $ is real-valued, 
or Fi(z) = (sin7)(l — cos7) _1 F 2 (z). In all cases, for z on the torus, 
Fi(z) = if and only if and F 2 (z) = 0, as required. □ 

5.6. Honeycomb frameworks and the kagome net. The explicit 
formulation of the matricial symbol function makes clear a simple ad- 
ditive principle which is useful for calculation : adding internal edges 
to a unit cell results in adding the same number of rows to the symbol 
function to form the new symbol function, and each of these rows has 
the standard form with constant entries. 

Consider the following honeycomb frameworks based on regular 
hexagons. 
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Figure 19. A 3- regular honeycomb framework. 




Figure 20. A 4- regular honeycomb framework. 



The 3-regular (symmetric) honeycomb framework has a 9 by 12 ma- 
tricial symbol function. There are vectors in the kernel of the multi- 
plication operator providing internal infinitesimal flexes and there are 
evident bounded infinitesimal flexes associated with "parallel deforma- 
tions" , for example. 

A 4-regular framework may be constructed by rigidifying the hexagon 
internal to the unit cell, by adding three cross edges. The matricial 
symbol function is correspondingly enlarged by three scalar rows, ap- 
pearing in rows 9 to 12 in the following matrix. 
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The cross-barred hexagon is continuously rigid but is infinitesimally 
flexible. We refer to this flex as the in-out flex. (See Whiteley [M] for 
related discussions.) It is straightforward to show that this flex extends 
periodically as a 1-cell-periodic infinitesimal flex of the framework. 

The motif matrix function is obtained by evaluating at z — w — 1. 
This has rank 8 reflecting the four independent 1-cell-periodic flexes 
corresponding to horizontal translation, to vertical translation, to local 
rotation of the cross-barred hexagon, and to the periodic "in-out" flex. 

We define the kagome net framework as the three dimensional crystal 
framework Qknet — {F v , F e ,7* 3 ) determined by the motif diagram in 
Figure 21. This figure shows a tetrahedron internal to a parallelepiped 
unit cell together with six additional non-internal edges, each of which 
extends a tetrahedron edge as indicated. 

w 

4 




Figure 21. A twelve-edged motif for the kagome net, Qknet- 



For the given motif the motif rigidity matrix R m {Qknet) takes the 
form of a 12 by 12 matrix valued function on the 3-torus. Different 
motifs, for the same translation group, give matrix functions that are 
related by permutations on columns and monomial multiplication on 
rows. It can be shown that the determinant of such a matrix function 
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is equal to a scalar and monomial multiple of the polynomial 

(z — 1) (w — 1) (u — 1) (z — w) (w — u) (u — z) . 

This is in complete analogy with the kagome framework and in the 
light of earlier discussions may draw the following conclusions. 

Theorem 5.13. (i) The kagome net Qknet is square- summably isostatic 
and possesses no internal infinitesimal flexes. 

(ii) The mode multiplicity function \x : T 3 — > Z + of the kagome net 
framework has support equal to the intersection ofT 3 with the six planes 

z = l,w = 1, u — l,z = w,w = u,u = z. 

5.7. Crystallography and Rigid Unit Modes. Rigidity theory and 
a Hilbert space operator viewpoint have led us to determine, ab initio, 
the matricial symbol function $>(z) of an abstract crystal framework 
and motif. The zeros of det <&(z) correspond to the phases of infini- 
tesimal periodic-modulo-phase wave flexes and the mode multiplicity 
function /x(u>) of $ detects the multiplicities of independent wave flexes 
for w. There are close connections between these mathematical obser- 
vations and the Rigid Unit Modes (RUMs) that are observed in certain 
material crystals through diffraction experiments. This connection is 
part of the motivation for the formulation of the function fi(z). It seems 
to us that the explicit algorithm for the passage from crystal motif to 
matrix function will provide a useful computational and theoretical tool 
for identifying RUMs and their relationships. 

Infinitesimal wave flexes for abstract crystal frameworks appear in 
the classical crystallography of Born von Karman theory. For exam- 
ple in the case of variants of quartz the low energy excitation modes 
result from rigid unit motions of SiO^ molecular tetrahedra within a 
tetrahedral net. These modes are therefore modelled by the first order 
periodic-modulo-phase infinitesimal rigid motions of abstract tetrahe- 
dral net frameworks. For this connection see, for example, Chapter 3 
of [35] . Comparisons of extensive simulations and experimental results 
have shown that the RUM modes of simulated crystals are closely cor- 
related with observed RUMs. For this see the seminal paper of Giddy 
et al [TT] . and also, for example, Dove et al [8] and Goodwin et al [H]. 

An interesting advance has been obtained recently by Wegner [32J 
who has derived mathematically, rather than through simulation, the 
RUM sets for various idealised tetrahedral crystals that model the ge- 
ometry of material crystals, including /3-cristobalite, HP tridymite, /3- 
quartz, a-cristobalite and a-quartz. These models are standard frame- 
work models that coincide with crystal bar-joint frameworks in our 
terminology and the results are obtained by determining vanishing de- 
terminants (zero sets of det $ in our formalism) by means of computer 
algebra and symmetry reductions. In particular Wegner obtains an 
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analytic derivation of some of the curious surfaces observed in experi- 
ments of Dove et al. 

REFERENCES 

[1] L. Asimow and B. Roth, The rigidity of graphs, Trans. Amer. Math. Soc, 245 
(1978), 279-289. 

[2] C.S. Borcea and I. Strcinu, Periodic frameworks and flexibility, Proc. R. Soc. 

A, doi:10.1098/rspa.2009.0676. 
[3] A. Bottcher, B. Silbermann, Analysis of Toeplitz operators, Springer- Verlag, 

Berlin, 1990. 

[4] R.Connelly, P.W.Fowler, S.D. Guest, B.Schulze, W.J.Whiteley, "When is a pin- 
jointed framework isostatic?" , International Journal of Solids and Structures, 
46 (2009), 762 - 773. 

[5] D. Cox, J. Little, D. O'Shca, Ideals, Varieties and Algorithms, Springer- Verlag, 
1992. 

[6] M. V. Chubynsky and M. F. Thorpe, Algorithms for three-dimensional rigidity 
analysis and a first-order percolation transition, Physical Review E 76, 2007. 

[7] A. Donev and S. Torquato, Energy-efficient actuation in infinite lattice struc- 
tures, J. Mech Phys. Solids, 51 (2003) 1459-1475. 

[8] M. T. Dove et al, Floppy Modes in Crystalline and Amorphous Silicates; in 
Rigidity Theory and Applications, Edited by M. F. Thorpe and P. M. Duxbury, 
Kluwer Academic/Plenum Publishers, New York (1999). 

[9] P.W. Fowler and S.D. Guest, A symmetry extension of Maxwell's rule for 
rigidity of frames, International Journal of Solids and Structures 37 (2000) 
1793 - 1804. 

[10] M. G. Frost and C. Storey, Equivalence of a matrix over 72. [s, z) with its Smith 

form, Int. J. Control, 28 (1978), 665-671. 
[11] A. P. Giddy, M.T. Dove, G.S. Pawlcy, V. Heine, The determination of rigid unit 

modes as potential soft modes for displacive phase transitions in framework 

crystal structures. Acta Crystallogr., A49 (1993), 697 - 703. 
[12] H. Gluck, Almost all simply connected closed surfaces are rigid, in Geometric 

Topology, Lecture Notes in Math., no. 438, Springer- Verlag, Berlin, 1975, pp. 

225-239. 

[13] I. Gohbcrg, P. Lancaster, L. Rodman, Matrix Polynomials, Classics in Applied 

Mathematics, SIAM, 2009. 
[14] A. L. Goodwin, M. Callcja, M. J. Conterio, M. T. Dove, J. S. O. Evans, 

D. A. Keen, L. Peters, M. G. Tucker, Colossal Positive and Negative Thermal 

Expansion in the Framework Material Ag3[Co(CN)6], Science 8 February 2008: 

Vol. 319. no. 5864, pp. 794 - 797. 
[15] J. Graver, B. Servatius and H. Servatius, Combinatorial rigidity, Graduate 

Texts in Mathematics, vol 2, Amer. Math. Soc, 1993. 
[16] S. D. Guest and J. W. Hutchinson, On the determinacy of repetitive structures, 

Journal of the Mechanics and Physics of Solids, 51 (2003), 383 - 391. 
[17] B. Grunbaum, G.C. Shephard, Tilings and patterns, W.H. Freeman, New York, 

1987. 

[18] K.D. Hammonds, M.T. Dove, A.P. Giddy, V. Heine, B. Winkler, Rigid-unit 
phonon modes and structural phase transitions in framework silicates, Ameri- 
can Mineralogist, 81 (1996), 1057-1079. 

[19] R.G. Hutchinson and N.A. Fleck, The structural performance of the periodic 
truss, Journal of the Mechanics and Physics of Solids, 54 (2006) 756-782. 

[20] A. B. Kempe, On a general method of describing plane curves of the nth degree 
by linkwork, Proc. London Math. Soc, 7 (1876), 213-216. 



46 



JOHN OWEN AND STEPHEN POWER 



[21] G. Laman, On graphs and the rigidity of plane skeletal structures, J. Engi- 
neering Mathematics, 4 (1970), 331-340. 
S. Lang, Algebra, Addison Wesley, 1993. 

J.C. Maxwell, On the calculation of the equilibrium and stiffness of frames. 
Philosophical Magazine 27, (1864), 294 - 299. 

J.C. Owen and S.C. Power, The non-solvability by radicals of generic 3- 
connectcd planar Laman graphs, Trans. Amer. Math. Soc, 359 (2007), 2269- 
2303. 

J.C. Owen and S.C. Power, Infinite bar-joint frameworks, Proceedings of the 
Symposium in Applied Computing, (SAC 2009) March 8-12, 2009, Honolulu, 
Hawaii, U.S.A. 

J.C. Owen and S.C. Power, Continuous curves from infinite Kempe linkages, 
Bull. London Math. Soc, 2009; doi: 10.1112/blms/bdp087 
J.C. Owen and S.C. Power, Frameworks, symmetry and rigidity, to appear in 
the Inter. J. of Computational Geometry and Applications, arXiv:0812.3785. 

E. Ross, private communication, July 2009. 

B. Roth, Rigid and flexible frameworks, American Math. Monthly, Vol. 88 
(1981), 6-21. 

W. Rudin, Functional Analysis, McGraw and Hill, NY, 1991. 
M. F. Thorpe and P. M. Duxbury, Rigidity Theory and Applications, Edited 
by Kluwer Academic/Plenum Publishers, New York (1999). 

F. Wegner, Rigid-unit modes in tetrahedral crystals, Journal of Physics: Con- 
densed Matter, Volume 19, Issue 40, (2007). 

S.A. Wells, M. T. Dove, M. G. Tucker and K. Trachenko, Real-space rigid- 
unit-mode analysis of dynamic disorder in quartz, cristobalite and amorphous 
silica, J. Phys.: Condons. Matter 14 (2002), 46454657. 

W. Whiteley, Matroids and rigid structures, in "Matroid Applications" ed. N. 
White, Encyclodedia of Mathematics and its applications 40 (1992), 1-51. 
B. T.M.Willis and A.W. Pryor, Thermal vibrations in crystallography, Cam- 
bridge University Press, 1975. 

M. Wyart, Rigidity-based approach to the boson peak in amorphous solids: 
from sphere packing to amorphous silica, preprint 2008, arXiv:0806.4596vl. 

D-Cubed, Siemens PLM Software, Park House,, Castle Park, Cam- 
bridge UK 

E-mail address: owen.john.ext@siemens.com 



Dept. Math. Stats., Lancaster University, Lancaster LAI 4YF, U.K. 
E-mail address: s.power@lancaster.ac.uk 



